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Abstract. Nonequilibrium collective motion is ubiquitous in nature and often results in a rich 
collection of intringuing phenomena, such as the formation of shocks or patterns, subdiffusive 
kinetics, traffic jams, and nonequilibrium phase transitions. These stochastic many-body fea- 
tures characterize transport processes in biology, soft condensed matter and, possibly, also in 
nanoscience. Inspired by these applications, a wide class of lattice-gas models has recently been 
considered. Building on the celebrated totally asymmetric simple exclusion process (TASEP) and 
a generalization accounting for the exchanges with a reservoir, we discuss the qualitative and 
quantitative nonequilibrium properties of these model systems. We specifically analyze the case 
of a dimeric lattice gas, the transport in the presence of pointwise disorder and along coupled 
tracks. 

1. Introduction. The modeling of collective transport phenomena is an important chal- 
lenge in theoretical physics, with possible technological and interdisciplinary implications 
(e.g. in biology, soft matter and nanotechnology) . In this context, motion of molecu- 
lar motors [How], single-file diffusion in colloid supsensions in confined geometry [Lin], 
spintronic devices [Zut], have motivated the study of a class of simple one-dimensional 
stochastic models. Remarkably, and despite their apparent simplicity, these models have 
been shown to be characterized by rich physical properties like nonequilibrium phase 
transitions, mixed phases, traffic jams, which are of interest on their own right. Hence, 
the goal of this article is to convey a pedagogical review of these lattice-gas models and 
of the techniques which allow to unravel their intriguing features. 

For the case of biological engines, it is only recently that modern experimental tech- 
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niques [Meh] have revealed the detailed causes of sub-cellular motion and transport. 
Nowadays, we know that every use of our muscles is the collective effort of a class of pro- 
teins called myosin that "walk" on actin filaments. All proteins that convert the chemical 
energy of ATP (adenosine-triphospate) in a hydrolysis reaction into mechanical work are 
referred to as molecular motors. These biological "engines" are highly specialized in their 
tasks: ribosomes move along mRNA strands while translating the codons into proteins, 
dynein is responsible for cilia motion and axonal transport, and kinesins play a key role 
in cytoskeletal traffic and spindle formation [How]. Rather than on the exact details of 
the molecular structure and function of motor proteins (see e.g. [Schl]), we are interested 
in phenomena arising out of the collective interaction of many motors. Early research 
along this line was motivated by mRNA translation that is managed by ribosomes, which 
bind to the mRNA strand with one subunit and step forward codon by codon. To in- 
crease the protein synthesis many ribosomes can be bound to the same mRNA strand 
simultaneously. This fact might induce collective properties, the slow down of an incom- 
ing ribosomes (and of the protein synthesis) due to steric hindrance caused by another 
ribosome in front, as was first realized by MacDonald [Mac], who set up a theoretical 
model for the translation of highly expressed mRNA. 

Much in the spirit of [Mac] our theoretical analysis focuses on simple, yet essential, 
collective processes underlying the motion and neglect the chemical or mechanical details 
on the molecular level of motor steps. Within this approach, even simple stochastic models 
are found to exhibit rich collective behavior. 

We will start this article by revisiting (Section 2) the main properties of the celebrated 
totally asymmetric simple exclusion process (TASEP). In Section 3, we review the main 
properties of a recently introduced model where the TASEP is coupled to a process 
accounting for the exchanges with a reservoir. In the following sections (Sec. 4-7), we 
outline recent results on variants of such a model and therefore discuss its robustness. 
In particular, the case of a lattice gas of extended particles is considered in Sec. 4, while 
the effects of pointwise disorder is addressed in Sec. 5. Sections 6 and 7 are devoted to 
the transport properties along two types of coupled tracks. In Section 8, we conclude this 
work by providing an outlook and our conclusions. 

2. The TASEP. Biological engines, like ribosomes and kinesins, often move along one- 
dimensional tracks (e.g. mRNA strands or microtubules). Each individual "step" is the 
result of a fascinating sequence of molecular events and is an interesting research topic 
on its own right [How], but is not our scope here. In fact, our approach is to devise simple 
and analytically amenable stochastic models. By identifying and mimicking the main 
ingredients underlying the collective transport, our goal is to get collective phenomena 
which we expect to be (qualitatively) comparable to those observed. 

A first step along this line of research was done in a pioneering work by MacDonald 
et al. [Mac], where the the total asymmetric exclusion process (TASEP) was introduced. 
This has since been widely recognized as a paradigmatic model in nonequilibrium statis- 
tical mechanics. It consists of a one-dimensional lattice (Fig. [1] left) with N sites labeled 
by i = 1, • • • ,N and with a spacing of a = L/N. Here, L is the total length of the lattice 
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Figure 1: (Color online). (Left) Schematic model of TASEP (particles are injected with 
rate a, move exclusively to the right being subject to hard-core exclusion, and are removed 
with rate (Right) Phase diagram of TASEP in a, (3 phase space. One distinguishes 
the low density (LD, green), the high density (HD, blue) and the maximal current (MC, 
yellow) phases, as well as the first (a = /3 < 1/2) and second order transition lines (in 
bold). 

and, for technical convenience, is often set to 1 with the lattice spacing then referred to 
as e = 1/A^. Particles (mimicking the motors) have an extension of the size of the lattice 
spacing and are subjected to hard core exclusion due to steric hindrance. Along the track, 
particles can jump to their right neighboring site with rate r (set to 1 in the following) 
provided the latter is empty, whence the occupation number Ui of site i takes either the 
value or 1. In addition, one has to specify the dynamic rules at the boundaries of the 
lattice. Contrary to equilibrium systems, the latter have a crucial influence and are re- 
sponsible for boundary-induced phase transitions [Kru] . While one can consider periodic 
or reflective boundaries, the richest behavior for the TASEP is obtained in the presence 
of open boundaries: at the left boundary {i = 1) particles attempt to enter the system 
with a rate a, while they exit the lattice at the right boundary {i = N) with rate (3. 
This is equivalent to two additional sites i = and i = N + 1 at the boundaries, which 
are connected to the system by the bulk dynamics described above, and are constantly 
set to the density a and 1 — f3 respectively. Because it does not fulfill the detailed bal- 
ance, the TASEP evolves into a steady state where a non- vanishing current is maintained 
between boundaries. Upon varying these boundary conditions, TASEP was found to ex- 
hibit phase transitions which - following general theorems [Mer] - are not even allowed 
for one-dimensional equilibrium systems in the absence of long-range interactions. 

2. 1 Analysis of the TASEP. While a severe difficulty for the analysis of non-equilibrium 
systems stems from the lack of universal concepts like the Boltzmann-Gibbs ensemble, 
models like the TASEP and its generalization can still be studied in many details by 
means of methods which are discussed in the remainder of this review article. Here, for 
the sake of clarity, we outline the analysis of the TASEP and some of its main statistical 
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properties. The latter are related to the density and current distributions, which are 
quantities of direct biological relevance. 

At any given moment, the system can be found in a certain configuration C made 
up of the occupation numbers at each lattice site. After a stochastic event occurred, i.e. 
the jump of one particle to a neighboring site, the new configuration is C . By its very 
definition, the TASEP is a Markov process. Hence, the probability of finding the system 
in the configuration C at time t obeys the the master equation: 

^P{C; = E i^'c'^cPiC; t) ~ ujc^c'PiC; t)] , (1) 

C'^C 

where the toc^c is the transition rate from the configuration C to C . Using, for instance, 
the so-called "quantum Hamiltonian formalism" [Sch], one can derive from ([1]) the equa- 
tions of motion for the density and current of particles. In fact the density of particle at 
site i, (rii), is found to obey 

^{n,) = - rii)) - - rii+i)), (i = 2, . . . , iV - 1), 

dt 

^{n,) = a[l-(ni)]-(ni(l-n2)), (2) 

37(n-7V> = {uN-iil - tin)} - P{nN), 
dt 

where the brackets stand for the average over the histories and, for brevity, we have 
omitted the time-dependence. From the discrete continuity equation (with Vj; = ji — 
ji-i), -^{ni) + Vj'i = 0, one also infers the expression for the local current of particles: 

ji = {ni{l - nj+i)) (3) 

The interpretation of (|2l3p is straightforward, but the analysis is difficult because these 
equations underlie an infinite hierarchy of equations: the density is connected to the two- 
point correlations, which depend on the three-point and so on. Despite these difficulties, 
an exact solution has been achieved by using the so-called matrix and Bethe Ansatz 
methods [Deri, Der2, Sch]. As these exact methods are mostly tailored for the TASEP, it is 
insightful to outline an alternative approach. The latter relies on a mean-field assumption 
and, while it is approximate in its essence, it has a broad range of application and allows 
to derive the main properties of the models discussed in this review. 

We now focus on the stationary properties of the TASEP and show how the phase 
diagram can be obtained within a simple mean-field theory. The latter relies on neglecting 
all the spatial correlations: {fnfij) = QiQj, where we have introduced Qi = (rii). In the 
steady state (dgi (t)/dt = 0) , this approximation leads to the following equations of motion 
for the density in the bulk: 

= ft_i(l - ft) - ft(l - ft+i), (4) 

while the expression for the current reads 

j^ ft(l - ft+i). (5) 

Clearly, the hierarchy of equations is now closed and, together with boundary conditions, 
forms a nonlinear set of equation giving rise to nontrivial solutions. A crucial observation 
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is the local conservation of the stationary current through the bulk. In fact, from Eq. (j4|) 
it follows that ji = ji-i- Actually, one remarks that such a relation is also embodied in 
Eq. ([2]) and holds true beyond the mean- field approximation. At this point, it is useful 
to adopt a continuum limit which turns the spatial lattice variable continuous. We thus 
consider a large number N of lattice sites and rescale the length of the system to L = 1, 
with infinitesimal lattice space e = L/N = 1/N <C 1. Introducing the spatial variable 
X — i/N, with < a; < 1, one can proceed with a Taylor expansion of the local density: 
Q{x±e) — Q{x)±edxQ{x) + ^e'^dlQ[x) + 0{e^). In the continuum limit and by neglecting 
the higher order terms in e, Eq. Q is recast as a simple first-order differential equation: 

{2q - l)9,g - . (6) 

The latter has to be supplemented with the boundary conditions gi(0) = a and q{1) = 
1— /3. Because we have a first order differential equation that needs to satisfy two boundary 
conditions, we are evidently concerned with an over-determined boundary value problem. 
There are three solutions to ([6]): Qh\i\k{x) = 1/2 does not satisfy either boundary condition 
(except for the special case a = (3 = 1/2), while q{x) = C can satisfy either the left or 
the right boundary condition, resulting in Qa{x) = a and Qf3{x) = 1-/3, respectively. To 
obtain a solution satisfying the boundary conditions, both solutions need to be matched. 
To resolve this issue, both solutions have to be considered valid in non-overlapping areas 
of the system. At the point where the current imposed by the left and right boundaries 
coincides (conservation of the current), those areas border and are connected by a sharp 
domain wall, such that g — Qa for < a; < and g ~ Qp for < < 1 ■ The domain 
wall is found to move with a velocity V — (3 ~ a [KoU] These simple considerations, 
together with the conservation of the current, have important implications and result in 
three distinct phases : 

• For (3 < 1/2 and a > (3, the current ja = Qa{l — Qa) of the left density solution 
Qa is higher than the current in the right part of the system and the particles are 
transported faster to the DW from the left end, then they head on the right. Hence, 
the domain wall is shifted to the left and the whole bulk density takes the value 
of Qj3 (except for a small boundary layer, which vanishes in the limit N — )■ oo). 
This situation therefore corresponds to a high density phase (HD), where both the 
density and the current are determined by the exit rate (3: the mean-field theory 
gives a spatially constant density pi — \ ^ (3 larger than 1/2 and a constant current 
ji — (3{1 — [3). In this case, the current is dominated by the low exit rate which acts 
as a bottleneck for the transport. 

• Along the same lines, for a < 1/2 and a < [3 the low entrance rate is the limiting 
factor for the particle current which is now given by ji = a{l — a). Since pi ~ a in 
always smaller than 1/2 in this parameter range, the phase is also termed the low 
density (LD) phase. 

• If both a and (3 become larger than the critical value 1/2, the current saturates 
and the density becomes constant p* = 1/2 independently of the parameters at 
the boundaries. The current is limited by the particle exclusion in the bulk and its 
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maximal value is j* — 1/4; therefore, this phase was termed maximal current (MC) 
phase. 

The (first-order) transition line a — (] < 1/2 separating the LD and HD phases 
requires special treatment. There, the velocity of the domain wall is zero, but Monte 
Carlo simulations and rigorous results [And] show that the DW actually performs a 
random walk in a domain with reflecting boundaries and an average over a sufficiently 
long time will therefore result in a homogeneous probability density over space (instead 
of a localized shock) g{x) = a + {1 — P — a)x. 

The above discussion is summarized in the phase diagram of the TASEP (Fig. [T]), 
where the LD and HD phases are separated by a first-order transition line, while the 
HD/MC and LD/MC phases are separated by a second-order transition hne. It has to 
be noticed that the phase diagram predicted by the mean-field theory coincide with 
that obtained by exact methods. Actually, the only features not accurately accounted 
for by the mean-field theory are the boundary layers and the density profile along the 
line a = /3 < 1/2 [Derl,Der2,Sch]. The somewhat suprising validity of the mean-field 
approach can be traced back to the density relation, j = g(l — g), which has been shown 
to be an exact result. 

3. The TASEP with on-off kinetics. With the aim to cover more realistic situations, 
there has recently been an upsurge of lattice gas models in the spirit of the TASEP. In 
fact, to devise more realistic transport models, some of the simplistic assumptions under- 
lying the TASEP had to be questioned. As an example, experimental observations have 
shown that ribosomes typically cover an area on the mRNA that exceeds the unit lattice 
space. To account for this situation the particles in TASEP have to extend over several 
lattice sites. This feature was accounted for by considering the £- TASEP, where extended 
particles (of length £) move unidirectionally being subject of excluded volume, and it was 
found that the TASEP phase diagram was not affected qualitatively [Sha]. It was also 
shown that in the presence of clusters of bottlenecks, mimicking slow codons along mes- 
senger RNA, may either give rise to local density perturbations or affect macroscopically 
the density profile [Chou,Tri]. 

Intracellular transport along cytoskeletal filaments has also served as a source of 
inspiration for driven lattice gas models. While in the TASEP model motors can only 
bind and unbind on the left and the right boundary respectively, cytoskeletal motors are 
known to detach from the track to the cytoplasm [How], where they perform Brownian 
motion, and subsequently reattach to the track. The interplay between diffusion in the 
cytoplasm and directed motion along the filament was studied [Lip] both in open and 
closed compartments, focusing on anomalous drift and diffusion behavior, as well as on 
maximal current and traffic jams as function of the motor density. In Ref. [Pari] it has 
been realized that, for an appropriate scaling of the on-off rates that ensures that particles 
always travel a finite fraction on the lattice, the on-off kinetics may not only give rise to 
quantitative changes in the transport efRciency but also to a novel class of driven lattice 
gas models. Actually, it was shown that the interplay between bulk on-off kinetics and 
driven transport results in a steady state exhibiting phase separation. In the following 



GENERIC PRINCIPLES OF ACTIVE TRANSPORT 



7 




^ (J mo • (J 



A 




Figure 2: Schematic model of TASEP/LK: the TASEP is extended by possible particle 
attachment and detachment in the bulk with rate cjajWd 

subsection, we review the main results of the studies [Parl,Par2]. 

3.1 The TASEP/LK. The crucial feature of cytoskeletal transport is the possibil- 
ity of bulk attachment and detachment with a finite residence time on the lattice. The 
latter can be understood as an effect of thermal fluctuations that may overcome the 
binding energy of the motors that is only of the order of several fceT. In this respect, 
attachment and detachment of motors is a stochastic process whose dynamic rules have 
to be defined. Parmeggiani et al. [Pari] chose to use Langmuir kinetics (LK), known as 
adsorption-desorption kinetics of particles on one- or two-dimensional lattices coupled to a 
bulk reservoir [Vill]. Particles can adsorb at empty sites and desorb from occupied sites, 
with microscopic reversibility demanding that the kinetic rates obey detailed balance 
leading to an evolution towards an equilibrium steady state (in the statistical mechanics 
sense). In this framework, the choice of LK is especially appealing as one is then faced 
with the competition of a prototypic equilibrium (Langmuir kinetics) and a paradigmatic 
non-equilibrium (TASEP) system. The resulting system - in the following referred to 
as TASEP/LK - is hence defined as follows: the well-known TASEP is extended with 
the possibility of particles to attach to the filament with rate lua and to detach from an 
occupied lattice site to the reservoir with rate wu- According to the type of ensemble 
(canonical or grand canonical) the reservoir is either finite or infinite. Here, the reservoir 
is assumed to be infinite and homogeneous throughout space and time. For a lattice ex- 
periencing only LK, the equilibrium coverage depends solely on the ratio K — lod/loa, 
the resulting steady state, given by Langmuir isotherm gi^ ~ K/{K +1), is completely 
uncorrelated for neglecting any particle interaction (except excluded volume). This ap- 
pears sensible as the diffusion in the cytoplasm is fast enough to fiatten any deviations 
from the homogeneous reservoir density. If we now consider the combination of TASEP 
and LK into the model displayed in Fig. [21 due attention has to be paid to the different 
statistical nature of the processes. On one hand, TASEP evolves into a non-equilibrium 
state carrying a finite current and, the number of particles being conserved through the 
bulk, the system is very sensitive to the boundary conditions. On the other hand, LK as 
an equilibrium process is expected to be robust against any boundary effects (especially 
for large systems). Combining both processes on the same time scale would unavoidably 
lead to a trivial domination of LK, as the bulk rates uja and ujd apply to a large number 
of bulk sites and would become predominant over the local entrance and exit rates a and 
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(3. To observe any interesting behavior (i.e. real interplay) between the two dynamics, 
one needs to ensure competition in the system. A prerequisite for the two processes to 
compete is to deal with comparable reaction rates, which requires an appropriate scaling. 
To this end TV- independent global detachment and attachment rates, flu and J^Aj are 
introduced, while the local on-off rates per site scale as 

= = ^. (7) 

The physical meaning of the above scaling is revealed by focusing on the involved time 
scales: a particle on the lattice performs a move on an time scale which is the inverse of 
the move's rate. Therefore, a particle spends an average time I/lob on the lattice before 
detaching. Bearing in mind that the TASEP jump rate is set to unity, a particle jumps 
to its right neighboring site typically after one unit time step. Therefore, the particle will 
travel a number A't = I/ij^d of sites before leaving the lattice. Compared to the lattice 
length this corresponds to a fraction of nx = Nt/N = l/(iVti;D)- To keep this fraction 
finite, and ensure effective competition between transport and on-off dynamics, t^D needs 
to scale as defined in Concretely, for the kinesin the ratio between (un)binding and 
diffusion rates is typically 10"-^ — 50~^ [How]. 

3.2. Mean Field Solution of TASEP/LK. To obtain density and current distributions 
of the TASEP/LK, we use again a mean-field approach. Hence, by neglecting any spatial 
correlations, along the same lines of Sec. 2.1, the stationary density is found to obey 

= - Qi) - - Qi+i) -ujDQi +wa(1 - Qi), (8) 

which has to be supplemented with the same boundary conditions as above. It is again 
useful to consider the continuum limit, where (at first order in e) the above equation 
reduces to 

i2g-l)d,g-LUT^g + LUA{l-g)^0, (9) 

with the boundary conditions g{0) = a and g{l) = 1 — /3. To proceed with the analysis of 
Eq. ([9]), it is useful to introduce the binding constant K = lua/ujb between the attachment 
and detachment rates. However, for the sake of clarity, here we restrict ourselves to 
the simple situation with same attachment and detachment rate, i.e. K = 1, and refer 
the readers to Ref. [Par2] for a comprehensive discussion of the mathematically more 
involved case where K ^ I. For equal on and off rates, ^Ia = ~ it is enlightening 
to rewrite Eq. ^ as {dxg — fl){2g — 1) = 0. Obviously, there are two solutions to this 
nonlinear differential equation: the homogeneous density gL — 1/2, given by the Langmuir 
isotherm, and the linear slope g{x) — Qx + C. The constant C being determined by the 
boundary conditions, which leads to one solution imposed by the entrance boundary, 
Qaix) = a + fix, and another one, namely gf3{x) = 1 — j3 — fl + fix, determined by the 
exit end. The complete density profile g{x) is the combination of one or several pieces 
of these three functions. Depending on how these are matched, we are led to distinguish 
several phases as explained below. 

3.3. Phase Diagram and Density Distributions. The only region in the parameter 
space of the TASEP/LK which is the same as for the simple TASEP is associated with 
fast entrance and exit rates. In this case, the system is again in a maximal current phase. 
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i.e. in a bulk-controlled regime, and the upper right quadrant of the phase diagram 
(Fig. [3l right) does not change compared to the phase diagram of the TASEP (Fig. [1]). 
In this regime, the additional bulk on-off dynamics at Langmuir isotherm gi^ = 1/2 does 
not result in any changes in the density distribution. In this case, non-equilibrium and 
equilibrium dynamics do not compete but cooperate. 

As in TASEP different solutions can be matched in various ways, the simplest being 
the connection by a domain wall between the left solution ga and the right solution g^. 
Depending on the current distribution, two possibilities have to be distinguished. As both 
solutions are spatially inhomogeneous, the corresponding currents ja and jp are strictly 
monotonic (Fig. [31 left). If the currents equal each other inside the system at a position Xw, 
the domain wall is localized at this position, as a displacement to either side would result 
in a current inequality driving the DW back to (see Fig. [3l left). As a consequence, 
the TASEP/LK displays multi-phase coexistence of low and high density regions (LD- 
HD phase) in the stationary state on any time scales. This has to be contrasted with the 
TASEP, where this behavior is only predicted for short observation times. Recently, this 
domain wall localization has been claimed to be observed experimentally [Nish]. 

If the matching of left and right currents is not possible inside the system, the known 
LD and HD phases are found. This is the case for one boundary condition being consid- 
erably larger than the other and depends on the slope of the density profile. This slope 
is determined by the ratio of the TASEP step rate and the bulk interchange rate fi: For 
large Q any density imposed by the boundary relaxes fast against the Langmuir isotherm 
of — 1/2, resulting in a steep slope of the density profile. 

This fact allows for the existence of two other phases with multi-regime coexistence. 
We could imagine a scenario in which the boundary imposed density solutions decay 
fast enough towards the isotherm to enable a three-regime coexistence of low density, 
maximal current and high density (LD-MC-HD phase). Furthermore, a combination of 
a MC phase with a boundary layer on one side and a LD or HD region of finite extend 
at the other boundary can be imagined. Not all these phases will be realized for every 
value of fi. Instead, the phase topology of two-dimensional cuts through the a, /?, fi-phase 
space changes. An example is shown in Fig. [3] (right). 

4. Dimers and Robustness. The results discussed in the previous sections have been 
obtained within the framework of generic models. The key notion is that a detailed 
description of a physical system is nonessential for the occurrence of the phenomenon, 
but a small number of ingredients encodes a mechanism that leads to the phenomenon. 
In order to corroborate that the models studied so far are representative for a wide class 
of models it is instructive to make a step towards reality and add new features to the 
model. 

In the case of molecular motors, their molecular structure already encodes how the 
motion of a single motors is performed. The most studied examples are kinesins, dyneins 
and myosin V which are known to consist of two sub-components referred to as lead and 
trail head. These heads bind specifically to a particular subunit of a long polar molecule 
acting as molecular track. As the dynamics of ribosomes along the mRNA motivated the 
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Figure 3: Reproduced from [Par2]. (Left) The DW connects the two densities pa and 
Pl3 {both dashed) and is localized at the point where the correspondent currents ja and 
match. Here low and high densities are connected together through a domain wall, and the 
resulting phase is termed LD-HD. Note the finite extend of the domain wall (localization 
length) that is only produced for Monte-Carlo simulations {solid wiggly line) and is not 
captured by mean-field results (solid line), reproduced from [Hinl]. (Right) Topological 
changes in the phase diagram of TASEP/LK for (a) n = 0.3, (b) n = 0.5, (c) n ^ 1. 




Figure 4: (Color online). Illustration of the model of driven dimer transport coupled 
to a reservoir 

study of ^— TASEP, for the transport of molecular motors, it was also suggested that a 
realistic minimal extension of the TASEP/LK should comprise dimers, i.e. objects that 
occupy two lattice sites at any instant of time. These extended objects move according 
the following elementary processes, depicted in Fig. |4l which underlie a generalization of 
the above TASEP/LK dynamics: 

Unidirectional hopping: If the lead head of a dimer occupies a site and the following is 
empty, the dimer advances one step to the right with unit rate. 

Entrance: At rate a, a dimer enters the lattice with its trail head provided that the first 
two sites are empty. 
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Exit: A dimer with its lead head on the last site leaves the lattice with rate (3 emptying 
the last two sites. 

Detachment: Everywhere in the bulk (i.e. trail head on a site i = 2, . . . , N — 2) a dimer 
leaves the lattice emptying two sites with a site independent detachment rate ud- 

Attachment: Everywhere in the bulk a dimer enters the lattice with its trail head, pro- 
vided that the considered site as well as its right neighbor is empty, with a site- 
independent attachment rate loa- 

Although at first glance it appears to be a straightforward extension of the monomer 
model, a theoretical description encounters a number of pitfalls. First, the mean-field 

approach based on ignoring any correlations between neighboring sites is too simplistic 
even on the level of pure on-off kinetics. There, the fact that two neighboring sites have 
to be empty to allow for an attachment process modifies the relation of the equilibrium 
density ph, i.e. the Langmuir isotherm, and the binding constant K = loaI^^d- Due 
to the hard-core exclusion, the corresponding statistical mechanics calculation reduces 
to a combinatorial problem [Tonk,Thom] yielding the thermodynamic equation of state 
Pj = PlW [McG]. 

The consequences of extended objects on the stationary density profiles for a 'pure 
TASEP have been studied in [Lak,Shaw]. In particular, they found that the phase diagram 
is in essence the one of monomeric TASEP, the differences being of quantitative nature 
only. A theoretical description has been achieved within a 'refined mean-field theory' 
that accounts approximately for the correlation induced by simultaneous hops of two 
monomers rigidly glued together. Special care had to be taken to link the density at the 
boundary to the entrance/exit rates. 

Recently, it has been recognized that the same refined mean-field ansatz is the proper 
starting point also to deal with the combined problem of driven transport (TASEP) and 
on-off kinetics [Pierl]. The Langmuir isotherm pL as well as an appropriate current- 
density relation is encoded in this approach. The density profile in the stationary state 
again results from a first order differential equation subject to two (non-trivial) boundary 
conditions. In particular, the scenario of phase separation and phase coexistence has been 
established for that model, and it has been rationalized that branch points and domain 
walls are generic features of such one-dimensional driven transport models. An important 
conclusion of these studies is that the phenomena that emerge upon weakly coupling a 
bulk process to a boundary controlled system are robust. 

Stochastic simulations have corroborated the analytically obtained phase diagram 
quantitatively. In particular, as a sensitive test of the theory the position of the domain 
wall has been measured and perfect agreement to the theoretical prediction has been 
found. Again it appears that a mean-field type approach is able to capture exactly the 
phase diagram, as well as the shape of the stationary density profiles. The weak coupling 
to the reservoir is sufficient to wash out the dynamic correlations that build up as dimers 
walk through the system. Similarly, the slow time scale for the on-off kinetics that intrin- 
sicly would lead to the random sequential adsorption is sped up by the continuous mass 
transport through the system. 
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5. Pointwise disorder. The presence of some form of randomness on the track is another 
reahstic ingredient which is highly desirable to incorporate in a model for intracellular 
transport. In fact, disorder may biologically be mediated by structural imperfections of 
the microtubular structure or proteins associated to the microtubules that change the 
affinity of the motors with the track. In the TASEP/LK language, imperfections on the 
microtubule can be mimicked by the presence of pointwise disorder [Chou] . In addition, 
studying the effect of disorder on dynamic and steady state properties of nonequilibrium 
systems is a major challenge in statistical physics [Stin]. 

In this framework, we have studied the influence of a bottleneck on the stationary 
transport properties of the TASEP/LK which is arguably its simplest, yet nontrivial, 
disordered version. Namely, as in the previous sections, we have considered the compe- 
tition between the totally asymmetric exclusion process and Langmuir kinetics in the 
mesoscopic limit [see Sec. 3.1, Eq. ([7|)] and in the presence of open boundaries, see Fig.[2l 
However, the system is now supplemented with the presence of a bottleneck of strength 
(7 < 1, which, at a given site k, slows down any incoming particles (Fig. [5l top). Depend- 
ing on the current of particles through the system and the strength of the defect, it is 
clear that the latter may have either a local (spike) or macroscopic effect (jump in the 
density profile) [Pier2]. 

In the region of the parameter space where the bottleneck is macroscopically relevant, 
the properties of the system are investigated through an effective mean-field theory in 
the continuum limit (k/N Xd). As illustrated in Fig. [H this approach is built on 
splitting the lattice into two subsystems on the left and right of the defect (respectively 
termed L and R). To ensure the local conservation of the current (which holds also in 
the presence of a bottleneck), one has to couple the L and R subsystems through the 
effective entrance and exit rates acff = q/{\-\-q) and j3cG = l/{l-\-q) on the right and left 
ends of L and R (Fig. [5]). Therefore, one is left with two coupled TASEP/LK processes 
whose current profiles have still to be matched. As the current has a spatial dependence 
(due to the Langmuir kinetics), the defect depletes its profile within a distance ^, that 
we call screening length. As the latter quantity increases with the strength of the defect 
and decreases with the attachment-detachment rates, it is an important concept for our 
analysis. In fact, depending on the screening length and the position of the defect, one 
has to distinguish between various scenarios, which are efficiently characterized by the 
corresponding carrying capacity. This notion stands for the maximal local current that 
can flow through the bulk of the entire system (including the defect). 

For the sake of simplicity, here we focus on the situation with equal attachment and 
detachment rates (i.e. 0,o = = As the screening length either covers entirely 
the system or part of it, we identify four types of carrying capacities Ci(x), z = 1, . . . , 4. 
Using the concepts of screening length and carrying capacity one readily obtains the 
phase diagram of the TASEP/LK system perturbed by a localized bottleneck. In fact, 
when the boundary currents dominate, the phase behavior of the defect-free system is 
recovered. Also, above some critical entrance and exit rates, the system transports the 
maximal current (carrying-capacity of the system), independent of the boundary rates. 
Between these two extreme situations, we find several coexistence phases: above some 
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Figure 5: Reproduced from [Picr2]. Schematic illustration of the effective mean- field 
theory: At the defect's location, i — k, the lattice is split into two subsystems (L and 
R), which are connected through the effective rates acs and /Jeff - The last step represents 
the continuum limit. The density p{xd) at the bottleneck's location is 1/(1 — q) on L and 
q/{l — q) on R. 



specific parameter values the system is dominated by the defect and the phase diagram 
is characterized by new bottleneck phases (BP). When the screening length spans the 
whole system, the associated carrying capacity, say Ci, never reaches the maximal value 
i* — 1/4: and one finds four new mixed phases. Among the latter, three phases (termed 
LD-BP, BP-HD and LD-BP-HD, see [Pier2]) are characterized by localized shocks in the 
density profile. When the screening length is short and the related carrying capacity, 
say C2, saturates at value j* within both subsystems L and R, one identifies nine novel 
bottleneck-induced phases. When the screening length is asymmetric (carrying capacities 
C3 and C4), i.e. covers entirely the subsystem R (L) and partially L (R), the are six new 
mixed phases [Pier2]. 

Each of these new mixed phases have been quantitatively studied within our mean- 
field theory and the analytical results have been checked against numerical simulations. 
As an illustration, let us consider the phase corresponding to the density profile of Fig. [HI 
In this case, the screening length is short and the carrying capacity is of type C2. As 
shown in Fig. [6] (right), near the extremities of the track, the current is imposed by the 
boundaries (j — ja and j — j^) and is below the maximal value j*, which is reached 
within the bulk. This gives rise to a low and high density (sub-)phases. At sufficient 
distance from the ends and from the defect, the current saturates at its maximal value 
J*, which corresponds to a maximal current (sub-)phase. Within a distance ^ from the 
defect (screening length; Fig. [6l right), the latter dictates the current profile and the 
system is in a bottleneck phase. As a consequence, the density profile displayed in Fig. [6] 
(left) corresponds to a phase termed LD-MC-BP-MC-HD. 

While the simple TASEP is known to be not affected qualitatively (only some tran- 
sitions lines are shifted in the phase diagram) by the presence of a localized defect 
[Jan,Kol2], the effective competition between TASEP and LK turns out to be very sensi- 
tive to the presence of pointwise disorder. We expect the same holds for other disordered 
versions of the TASEP/LK, e.g. in the presence of clusters of bottlenecks [Chou], to which 
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Figure 6: Example of the TASEP/LK density (left) and current (right) profiles in the 
presence of a bottleneck located at Xd — 1/2. The system (with parameters q = 0.3, 
a = P = 0.2, VLd = = l-^ and N = 4096) is in a mixed phase (LD-MC-BP-MC- 
HD), where (from the left to the right end) a low density (LD) and maximal current 
(MC) subphases are followed by a bottleneck- induced (sub)phase (BP). The latter is 
also concatenated to a MC and a high density (HD) subphases (see text). Results of 
simulations (continuous line) are compared to analytic predictions (dashed line). 




Figure 7: Reproduced from [Rei]. Illustration of a two-lane exclusion process. Lane I and 
II possess individual entering rates, resp. a^^ as well as exiting rates, (3^ resp. (3^^ . 



the above effective theory can be adapted. 

6. Multiple parallel lanes. Microtubules, the intracellular tracks for molecular motors 
like dynein or kinesin [How], are built of 12-14 parallel lanes. Although it has been 
revealed that the motors proteins typically remain on one track while proceeding on the 
microtubule, the statistics of deviations (random lane changes) is so far unknown. In the 
following, we aim to investigate effects of a small, but finite number of lane changes on 
the resulting density profiles along the lanes by studying a transport model consisting of 
two parallel lanes I and II, see Fig. [71 

Models with indirect coupling of lanes or strong direct coupling through frequent 
lane changes have been investigated in [Popl,Pop2,Pro,Mit], while the focus of the work 
presented in the following is on weak coupling emerging from only few lane changes. 
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Particles are injected at the left boundary at rates (a^^) on lane I (II), provided the 
respective lattice site is empty. In the bulk, particles hop forward at a constant rate 
which we set to unity; in addition, they may change from lane I to II and back at a small 
rate lu (weak lane coupling) , always under the constraint of simple site exclusion. Having 
reached the right boundary, particle extraction occurs at rates f3^ (Z?^^)- 

Interestingly, the above introduced model allows for an interpretation in terms of spin 
transport. Indeed, on a given lattice site, a particle may be on lane I or II, which may 
be viewed as an internal state, e.g. the spin state of an electron. Simple site exclusion 
translates into Pauli's exclusion principle: spin-up and spin-down may simultaneously 
occupy the same site, though not two particles of the same spin state. The resulting 
dynamics, depicted in Fig.[8l may have possible application within the field of spintronics 
[Zut]. 
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Figure 8: (Color online) Reproduced from [Rci]. The same model as above, in the spin- 
transport interpretation. Particles respect Pauli's exclusion principle; spin flips occur at 
a (small) rate uj. 



As for the problem of transport with attachment and detachment of hard-core particles 
in the bulk [Sec. 3.1, Eq. we employ a mesoscopic scaling for the exchange rate oj by 
keeping the gross lane change rate, defined as = cuL, as constant for large system size L. 
The latter yields a measure for the total number of lane changes that a particle undergoes 
while traversing the system, which is thus kept fixed, in this way ensuring a competition 
of lane change events with the boundary processes (injection and extraction). 

This model exhibits a variety of different phases for the emerging non-equilibrium 
steady state, controllable by the injection and extraction rates as well as the lane change 
rate. A typical scenario is presented in Fig. [9l where the language of spin transport has 
been adopted: denotes the average density of spin- up (i.e. the concentration of particles 
on lane I); and stands for the density profile of spin-down (lane II) particles. While the 
density of spin-down remains at low values, the density of spin-up undergoes a transition 
at Xuj from a region of low to a region of high value with the emergence of a localized 
shock. In the spin transport interpretation, we encounter a spontaneous polarization: the 
densities of both spin states rather equal each other to the left of the shock, but largely 
differ to the right, yielding a polarization there. A key point for the analytical description 
is the continuity (though not smoothness at Xw) of the individual spin currents, shown 
in the inset. 

In the same spirit of the previous sections, we have taken advantage of a mean-field 
approach and a continuum limit to analytically describe the stationary state density pro- 
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Figure 9: Reproduced from [Rci]. Typical scenario for the non-equilibrium stationary 
state: a shock forms in the spin-up density profile, inducing a spontaneous polarization. 
The spin currents, shown in the inset, are both continuous. Solid lines correspond to 
the analytical solution, while dashed lines indicate results from stochastic simulations on 
lattices with L = 2000 (blue or dark grey) and L — 10000 (red or light gray). 

files. Surprisingly, comparison with stochastic simulations (finite size scaling) uncovers 
the apparent exactness of the analytical results in the limit of large system sizes. Fur- 
thermore, our analytical approach allows for the derivation of phase diagrams, leading to 
the full knowledge of the system's phase behavior. 

7. SEP and TASEP coupled on a ring. While not only in the TASEP but also in 
all extended models presented above, the strategy of a mean-field approach and the con- 
tinuum limit has proven succesful to obtain correct phase diagrams and density profiles 
in the limit of large system sizes, this cannot a priori be taken for granted. This was 
recently demonstrated in a study on the competitive effects of driven versus diffusive 
motion in lattice gases [Hin2] . As prominent examples serve the two paradigms of driven 
and diffusive transport in lattice gases: the purely driven TASEP, as introduced above, 
and the purely diffusive symmetric exclusion process (SEP). The analyzed model systems 
is a periodic setup (Fig. \TU[ left) consisting of a SEP and a TASEP part of equal size N 
with Np particles inside the system. The hopping rate of unity on the lower lane signifies 
a choice of time units and the hopping rate D on the upper lane fixes the ratio between 
driven and diffusive motion. 

In a first step, we proceed along the lines of analysis used above to derive a mean- field 
phase diagram. Here, the central idea is to decouple the two parts of the ring and con- 
sider them as separate lanes with effective entrance and exit rates. Known results for 
both processes can then be applied. While TASEP 's non-equilibrium phase transitions 
have been introduced in Sec. 2, SEP is known to exhibit a linear density profile con- 
necting the boundary values [Sch]. The periodicity of the system has two important 
consequences: current conservation and conservation of particle density Hp — Np/2N. 
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The current conservation allows one to establish the relations a = D6 and j — 1 — P 
between the boundary densities. The system's state can then be described by the two 
remaining control parameters D and Up. The non-equilibrium steady states have been 
classified and a phase diagram is derived as a function of those two parameters. 
In the case of TASEP /LK a mesoscopic scaling was used to ensure competition for arbi- 
trary large systems. To enlight the behavior of the ring system, the current is of crucial 
importance. While the lower lane carries a current = p(l — p) independent of system 
size, the upper lane's current is determined by the boundary condition's difference as 
= (7 — S)D/N reminding of "Pick's law". Obviously, this current decreases with the 
number of sites, leading to a "frozen" system due to an unphysical vanishing of current. 
However, a mesoscopic scaling d = D/N ensures a competition between both processes 
and a finite current as expected in the case of continuous diffusion in the limit N 00. 
Using a mean- field analysis and the conservation laws, phase transitions can be identified 
at characteristic changes in the TASEP density profile. As for the simple TASEP a LD 
and a HD phase are observed. However, a crucial difference is the existence of a localized 
domain wall that is observed in simulations and whose position Xuj 00/11 be derived as a 
function of the new control parameters as = (—3 + 3d + 4np)/(4(i — 2). There exists 
thus a LD-HD phase of finite extend in the ring system. This can be explained by particle 
conservation and is a consequence of the periodic system being canonical instead of the 
grand canonical TASEP. Additionally, the variance of the domain wall position is found 
to decay much faster with N^^ compared with a N~^^'^ decay in TASEP. The extend of 
the LD-HD phase is computed to (see Fig. [101 right) : 

3-4np 

d = tor Xw = U , 

d = iup - 1 for Xw = 1 . (10) 

By the introduction of a finite LD-HD phase the first-order phase transition between LD 
and HD phase is substituted by two second order transitions. Similar to TASEP a MC 
phase can be identified with phase boundaries: 

d = and d = . (11) 

16np - 4 8 - 16np ^ ' 

In a second step, we used Monte Carlo simulations to validate the derived phase 
diagram. It is found that strong deviations (symbols in Fig. [TO] right) occur. Thinking of 
the success of the mean-field method in similar systems, this comes as a surprise at first 
sight. Investigating its origin, a closer look at the density profile reveals that the upper 
part of the density profile is actually strictly convex instead of linear. This is a result 
of a strong time scale separation of the two subprocesses introduced by the mesoscopic 
scaling. While diffusive hopping events happen on a very short time scale, the entrance of 
a particle at the r.h.s. junction is a rather rare event. Between these entrance events the 
particle surplus at the r.h.s. relaxes diffusively to the left. As a result the time averaged 
density profile appears convex and "Pick's law" for the current is no longer valid. The 
current is effectively diminished since it is now bounded to its lowest value at the point 
of smallest local gradient. This explains the qualitative shift of the critical point in the 
phase diagram to higher d. 
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Figure 10: {left) (Color online). Schematic illustration of the ring system. The dynamics 
of the upper lane is governed by a SEP with rate D, and the dynamics of the lower lane by 
a TASEP with unity jump rate. The boundary densities are indicated at the junctions. 
(right) Reproduced from [Hin2]. Phase diagram obtained by MP (solid lines) exhibits 
four phases. Simulations [MC (triangles) and LD-HD (circles) phase boundaries] reveal 
a failure of the MP analysis. 



The observed behavior bears analogies to forest fire models [Dro], where out of two 
processes of different time scales the faster one strives towards equilibrium and the slower 
one occasionally forces the system to a non-equilibrium state. It furthermore reminds of 
the coexistence line in simple TASEP with open boundaries where a linear profile is 
obtained as the time average over a randomly walking domain wall. Hence, the study of 
the ring system highlights the importance of time scale separation and points out the 
limits of mean-field theory in lattice gas models. 

8. Outlook and conclusion. In this work, we have described how complex nonequi- 
librium transport phenomena, e.g. the motion of molecular motors along niicrotubuli or 
transport of spins, have inspired stochastic models displaying rich nonequilibrium collec- 
tive properties. Namely, in the same spirit as mRNA translation was modeled by means 
of the totally asymmetric simple exclusion process (TASEP) [Mac], intracellular transport 
processes have motivated the study of a broad class of lattice gas models. To account for 
exchanges with a reservoir (e.g. mimicking the binding/unbinding of motors from/to the 
cytosol) and the track, TASEP was supplemented with an on-off (Langmuir) kinetics. An 
effective competition between driven motion and Langmuir kinetics arises provided the 
related attachment and detachment rates scale with the inverse of the system size. In this 
case, a finite residence time on the lattice is ensured and results in the emergence of a 
new class of non-equilibrium features, like multi-phase coexistence, localized shocks and 
enriched phase diagrams. Along these lines, we have reported recent results for physically 
motivated variants of the above model. 

Specifically, we have discussed transport properties of a one-dimensional gas of dimers, 
which was motivated by the fact that molecular motors like kinesin have two heads. The 
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finite extension of the particles then has led to a breaking of particle/hole symmetry 

accompanied by a failure of a naive mean-field approach. A more sophisticated approach 
had to be developed that allows to calculate the exact phase diagram (essentially) an- 
alytically. It has been argued that the scenario of phase separation is generic for this 
kind of models. We have also considered the effects of pointwise disorder in its simplest 
form, i.e. in the presence of a localized bottleneck slowing down any incoming particle. By 
means of an effective mean-field theory and numerical simulations, we have shown that 
such a single defect is responsible for the emergence of numerous new mixed phases and 
can result in the simultaneous formation of two traffic jams. The presence of multiple 
parallel lanes on which motors proceed can lead to further interesting phenomena. We 
have studied a two-lane model with rare lane switch events of the particles and discussed 
the mapping to a possible spintronic device. Depending on boundary rates, the station- 
ary density profiles can exhibit jammed states appearing on one or both lanes. Finally, 
the study of competition between driven and diffusive motion on two lanes in a periodic 
setup has shown the limitations of mean-field theory in the presence of strong time-scale 
separation. 
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